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a Kubo formula. We show how to sum an infinite series of planar ladder 



scattering processes between hard particles in the plasma. The result agrees 

with the one obtained previously from a kinetical approach. 
ll.10.Wx, 12.38.Mh, 05.60.-k 



*wtbmare j @lg . elm . es 
^wtpvabamOlg . elm . es 



I. INTRODUCTION 

Color conductivity has turned out to be a significant quantity for it enters directly into 
the theory governing the dynamics of long wavelength excitations in the QCD plasma. 
This was not obvious a priori, because a hydrodynamical description of the plasma usually 
involves transport coefficients associated with relaxation processes of momentum degrees 
of freedom for which the collision frequency is u p ~ g 4 T log(l/g) [|TJ, while the collision 
frequency for color relaxation is v c ~ g 2 T log(l/g) 0. However, a major achievement in 
understanding the dynamics of very soft excitations in the QCD plasma, has been Bodeker's 
finding of an effective theory described by a Langevin equation, in which thermal noise 
is parametrized by the static color conductivity H,§|. Although he did not use explicitly 
the concept of color conductivity, the leading logarithmic order result was implicit in his 
computation. As an intermediate step towards Bodeker's theory, Arnold, Son and Yaffe || 
have proposed a Boltzmann equation from which color conductivity can be obtained. Along 
the way, some points of previous studies of color conductivity [|,|| have been clarified. 
Later, a rigorous derivation of the Boltzmann equation, starting from the quantum field 
equations, was given in Ref . [|7| . The Boltzmann equation leading to Bodeker's theory (and 
color conductivity) can also be obtained starting from classical transport theory || or from 
the kinetic formulation of the hard thermal loop (HTL) effective theory ||. Furthermore, 
a recent diagrammatic analysis by Bodeker [|H|] of the polarization tensor beyond the hard 
thermal loop approximation has also displayed the form of the Boltzmannn equation. Some 
calculations of the color conductivity at next-to-leading log order have appeared over the 



last year [11,12 



These studies of color conductivity at leading-log order have used a kinetic approach. 
Other transport coefficients in hot gauge theories have also been studied over the last 
years |l| JT3|Jl^ J^1 using this approach. But there is an alternative way to compute trans- 



port coefficients, namely, through the use of Kubo formulas. However, this approach poses a 
difficulty: a summation of an infinite series of diagrams is required. Only for a scalar theory 



has this summation been carried out explicitly fL5 |. The higher order contributions were 



shown to come from ladder diagrams. In addition, the equivalence under certain conditions 



of a kinetic equation to the field theoretical computation was established [ p~5] , [T6|] . For the 
case of gauge theories it has been suggested that this kind of diagrams contributes at lead- 
ing order ||. However, because of the difficulty in carrying out these computations, this 
approach has not been pursued further. 

The purpose of this paper is to show how the ladder diagrams can be explicitly summed 
at leading-log order within the imaginary time formalism of thermal field theory. We will 
work at high temperature so the coupling constant is small and all zero temperature masses 
can be neglected. 

II. COLOR CONDUCTIVITY 

The static color conductivity associated with color flow may be defined by the constitutive 
relation (j*) = a l J b El, where El are the components of an external, constant chromoelectric 
field and (j*) is the ensemble average of the spatial part of the conserved color current 
density. Using linear response theory one may express color conductivity in terms of the 
low frequency, zero momentum limit of the retarded correlation function between color 
currents |18|| . The corresponding formula is 

< = -^imng a V,q = o)| 9 o =0 . (i) 

Because of isotropy, the dependence on the spatial and color indexes is very simple, of? = 
Sij S ab a c . The most efficient way to compute a c from the above Kubo formula is to exploit 
the Lehmann representation, which provides a direct connection via analytic continuation 
between the retarded Green's function and its counterpart in the imaginary time formalism, 
U R (q°, 0) = U.(v q = ~iq° + 0+ 0). 

The dominant contribution would come from the diagrams in Fig. [I] with all other dia- 
grams of higher order in the loop expansion being subleading a priori. However all ladder 



diagrams of Fig. |2] will contribute at the same order, as will be shown in our explicit com- 
putation. This same kind of diagrams has been shown by Jeon |15|] to give contributions 
of the same order in A as the one-loop graph for the viscosity in A0 4 theory. In the con- 
text of gauge theories, it has been suggested || that the same kind of diagrams gives the 
leading order contribution. This can be understood in the following way. In the limit of 
zero spatial momentum and small frequencies, there are going to be pairs of propagators 
which will carry the same momentum, and this would lead to an on-shell singularity for each 
product of propagators with nearly the same momenta. The solution is to include in the 
propagators of the side rails their thermal widths I\ Thus, for each pair of propagators, the 
singularity is replaced with a term 1/T. Hence, for a general ladder graph with n rungs we 
will have g 2 (l/T) n+1 (g 2 ) n where the first g 2 comes from the external vertices and each rung 
introduces a factor g 2 [[H||17]]. For gauge theories, the thermal widths for hard momenta are 
proportional to g 2 Tlog(l/g), where the logarithm arises from the infrared behavior of the 
theory. As we shall see, a similar logarithm is produced by each integration over the small 
momentum transfer in the rung. So, both the one-loop and the ladder graphs are propor- 
tional to g 2 T 2 q°/T. Other ladder diagrams with two gluons lines in the rung are subleading 
because the four gluon vertex has one more power of the coupling constant than the three 
gluon vertex. 

Although we lack a rigorous power-counting scheme to dismiss all other possible higher 
order diagrams in the loop expansion, we do not expect them to contribute at leading order, 
since we are not aware of any other way in which the coupling constants of each new vertex 
could be compensated. 

We now turn to the main point in this work, that is, how to sum the ladder diagrams of 
Fig. 0. First, we study the case when the side rails of the ladder are quark propagators. 



A. The quark contribution 

To include the thermal width in the propagators of the side rails, we modify the quark 
propagator by changing the delta functions of the spectral density for the free propagator 
by Lorentzians of width 2 r )f 

p(u 7 p) = S+(P)- /*[( 2 + S.(p) .^L. 2 , (2) 



where «S±(p) = (7 =F 7 • p)/2 and 7/ is the fermion damping rate at hard momentum [19 



N 2 - 1 q 2 
If = — — TWl/o). (3) 
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From this spectral function the fermion propagator follows, as 

SVn,p) = r ^L£^L = 5 + (p)A + (p) + SL(p)A_(p) 

J -00 2/K %UJ n — OJ 



S+(P) + ^~(P) / 4 N 



iuo n - |p| + ry/-sgn(covJ icu n + |p| + Z7/Sgn(a;, 
To sum the ladder contributions we write the following equation for an effective vertex 

d 3 k 



r> a (p + q , P ) = r d a + g 2 TY,J 



(2tt) 3 
x 7 a t b S(k + p + q)r l > a (k + p + q,k + p)S(k + p^t'D^k), (5) 

where the zero components of the momenta are % times the corresponding Matsubara fre- 
quencies. To" = 7*t a is the tree level interaction vertex and t a are the SU(N C ) generators 
in the fundamental representation, with tr(t a t b ) = 5 ab /2. The gluon propagator D a p is the 
HTL resummed propagator whose inclusion will be shortly justified. This equation is shown 
diagramatically in Fig. |3|. To solve it exactly is beyond all hope, so we will try to extract 
the dominant behaviour by making some approximations. 

The inclusion of the HTL propagator in the rung is based on the fact that the interaction 
of two particles must be modified due to many body effects for small momentum transfer 
k <T. It has the form 

A* (*) = v U^ ul + ^ +YiT{k) + *^ | k | a + n £ (fc) ' (6) 
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where V^(k) = 5ij — kikj, Vj a = Pj 4 = 0. In order to perform the Matsubara sum in 
Eq. (fJD, it is convenient to introduce the spectral representation of this propagator 



n / IN l f f f 00 C?CU 1 

L> a/ j(u; n ,k) = T^d^d^ + / t^ - " 

K 2 J-oc ZTT %ijj n — UJ 



-p T (u, k)VL + Pl(u, k)&ta£ 



1/3 



(7) 



and a double spectral representation for the vertex, whose general form parametrised by two 
functions Fl a 2 (uJi, 0J2) is given by 



r> a ( P + q, P ) 



00 duJ\ du 2 
-00 2tt 2tt 



F{ a {oo u u 2 ) 
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F^(uj 1 ,UJ 2 ) 



(8) 



Inserting the spectral representations into Eq. (|5]), this reduces in the high temperature limit 
uj <^T to 

duj <i 3 k 



r l,a (o;p + u q , uj p ; p) = r ' a - g T 



- 1 a t b S{u p + v q ^ + V ) 



2n (2tt) 3 
x r> a (u p + u q , uj p - k + p)S(uj p , k + p) 7 /3 t b P ^ (a,,k) , (9) 

where the factor T/u; comes from the term containing the Bose-Einstein distribution N(uj) 
produced in the summation and uj p , u q denote the fermionic and bosonic Matsubara frequen- 
cies of the external lines. Thus, the summation has been traded for an integration over 

UJ. 

To proceed further in the integration, let us recall the relevant kinematical domain from 
which arises the logarithmic sensitivity of the hard damping rates to the magnetic mass of 
order g 2 T. As we shall see, the leading-log order in the color conductivity results from the 
momentum transfer in the same domain. 

As is well known, the transverse part of the interaction is dynamically screened due to 
Landau damping of the magnetic virtual gluons in the regime \u\ < |k|, and static Debye 
screening cuts off effectively the small angle divergences of the longitudinal interaction. In 
terms of the gluon spectral density, the (k, uS) integral expressing the hard thermal damping 
rate (|]) is ||l9H2lH given by 



N} 



If 
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2N C 4tt j Aii 



dk k I - — \p L (u, k)+{ 1 - -pr )pt{u, ■ 

-k 2lTUJ k 



(10) 



where the lower cutoff A m j n , which is taken of the same order as the hypotetical magnetic 
mass ~ g 2 T, regularizes the infrared divergence of the integrand. The leading-log contri- 
bution to the damping rate arises from the exchange of quasistatic magnetic gluons in the 
uj <C k limit. The reason for this is that the denominator of the integrand kpr/uJ vanishes 
in the static limit for k = 0, while the denominator of kpi/u remains finite in this limit. In 
that domain, we can use the approximate expressions for the spectral densities 

Mm 2 D uk 

PL ^ k) = k(k 2 +m ir (12) 

where m 2 D = (Nf + 2N c )g 2 T 2 /6 is the Debye mass. Inserting them into Eq. (pTOj ) produces a 

dominant logarithmic term 

N 2 -lg 2 T r^dk I ' nm 2 D \ N 2 - 1 g 2 T f A ^ dk 

If = — ^r / -r arctan -—%• w — c - — / — (13) 

u 2N C 2vr 2 A min k \Ak 2 ) 2N C 4vr A min k V ; 

coming enterely from the transverse part, together with a finite term coming from the 

longitudinal part which is of order g 2 T. Here, we only need the most singular logarithmic 

term, but for a computation of other transport coefficients, all contributions are possibly 

required. The upper cutoff A max must be chosen of order gT because when k 3> gT, the 

contribution from hard momentum transfer k ~ T is computed using a non resummed gluon 

propagator. Thus the contribution to the damping rate from this domain is proportional to 

9 * r3 f % (14) 



If the upper cutoff were chosen of order T in Eq. ([13|) , one would make the mistake of relying 
on an integrand which in the range gT < k < T does not contain the correct dependence 
on the momentum transfer. So the result would be wrong by a factor of two. (The correct 
spectral density in that range is proportional to m 2 D u>/(k 3 (k 2 — uj 2 )).) 

Briefly, concerning the (k, uj) integration in the vertex equation ([5|), this can be computed 
by retaining only the transverse part, giving the result 

rk ^PL ~ 1 (15 ) 



This result is the same we would have obtained if, in the initial vertex equation fl5|), we 
had retained only the static mode in the sum with the substitution D a p — ► (5ij — kikj)l/k 2 . 
The following integration in k must be cut off by a semihard scale A max ~ gT on one 
side and a soft scale A m j n ~ g 2 T on the other side. In addition, for a pair of fermionic 
external lines nearly on shell, for example p° ~ | p | , p° + q° ~ |p|, we have only to retain 
in the vertex equation @ one of the nearly singular products, A + A + in that case. Of 
course, there is another pair of values of the external lines leading to a nearly singular piece, 
p° p« — |p|,p + q° ~ — |p|, for which only the product A_A_ must be retained. This 
additional simplification is due to the fact that when k <C p, a small momentum transfer 
cannot change the mass shell of the hard particles involved in the scattering. 
Now, we make the following ansatz for the efective vertex [] 

P'°(p + q,q) = (AY + 5 7 V + C7W> a , (16) 

where A,B,C are some scalar functions which can depend on p°, q° and |p| . First we consider 
the term with A + A + . The equation for the efective vertex is now 

V>\p + q,q)= IT - ^- jj^y 3 A + (p + k + q)A + (p + k) 

x (A + B + C) ( 7 ° - 7 " k(pTk) • k) (pTfcj'jl. (17) 

|k| 

The next step is to compute the integrals over the momentum transfer k. It is useful to 

introduce a momentum variable u, defined as u 2 = (k + p) 2 , so that 

dkk I duu . . . (18) 



(2tt) 3 '" (2n) 2 p J A min ' J\k- P \ 

One may decouple the limits of integration in the following way. In the limit q = 0, q° — > 0, 
with the pair of external fermion lines nearly on-shell at hard momentum, p° ~ |p| ~ 
T, we are concerned with a pair of internal fermionic lines sharing nearly the same loop 
momenta, so we can anticipate that values of u ~ p° and u ~ p° + q° are important in the 



1 The posible structure of the vertex is greatly simplified by the fact that q = 0. 



integration process. Since most of the contribution is in the region u ~ p°, the limits on 
the integral over u can be extended to ±00, with the contribution outside the initial range 
being small. Furthermore, if the scalar functions A, B, C are not strongly dependent on p 
we can take them out of the integral, along with the factor u, evaluated in u = |p|. Given 
these approximations, the u integration results in 

r du l \ = 2nV(u Jp ,u q ) ^ 

J -00 itjp — u + ?7/sgn(cJ p ) iuj p + iu q -u + ?7/Sgn(u> p + v q ) \u q \ + 27/ 

where T>{uj p , u q ) = 9(—u p )9(lj p + u q ) + 9{ui p )6{—uj p — u q ). We see that the thermal width 

regulates a singular behaviour when v q vanishes. Upon completion of all these integrations, 

it is straightforward to solve for the scalar functions. Finally, we obtain for the vertex 

function when the external fermions are nearly on shell at leading logarithmic order 



2j f V(u p , V q 

(JV?-l)|i/,|+2JV? 7/ 



1 (lu p + v q , u p , p) = 7 t - j--z — - — „ — 7 V t , (2(Jj 



where the logarithmic dependence on g coming from the k integral has been replaced by the 
same logarithm from the fermion damping rate. This result is consistent with the assumption 
about the smooth dependence on |p| of the scalar functions of the vertex. If we consider 
the term A_A_ we arrive at the same result. Here it should be emphasized that the vertex 
correction arising from ladder summation is of the same order as the one at tree level. 
It remains to compute the polarization tensor 

nS(^O) =^TX)/^tr[P> I , + i/ s ,a; 1 ,;p)5(a; 1 ,,p)yt 6 5(a; p + i/ < ,,p)] . (21) 

The most difficult task is to do the Matsubara sum. It will be useful to introduce a double 
spectral representation of the product of Green's functions as follows. Let A(p + q,p) = 
G(p + q)G(p), where the Green's function G(p) has a single spectral representation 

= r» ^/^p)_ (22) 

KFJ J -00 2tt p° - u y ' 

Then, it is easy to check that A(p + q,p) admits the double spectral representation 

Fi(wi,w 2 ) F 2 (uj 1 ,uj2) 



1 dui duj 2 
A(p + q,p) - / 

-00 Z7T Z71 



(p°-^)(g -cu 2 ) {jfi + gO- Ul ^gO- U2 ) 



(23) 



where the corresponding F functions are 

F t (u 1 ,u 2 ) = p{u>x + u 2 ; p + q)p{ui, p), (24) 

F 2 (ui,u 2 ) = -p(uJi; p + q)p{uJi - uj 2 ; p). (25) 

Similarly, the product of the vertex correction and two propagators, 

SA(p + q , P ) = ^4-. -1-. --, : pJL ( 26 ) 

+ c\v q \ iuo p - |p| + 27 sgn^p) iuj p + «z/ g - |p| + 27 sgn(cu p + v q ) 

has the same representation with the 5F functions 

SF( x = 2fr(7 2 + |p| 2 + uf- 2|p| coi - |p| uj 2 + iOjUh) - c^u\ 

U U 2) (7 2 + K-|p|) 2 )(7 2 + (a;l-|p|) 2 -2|p|^ 2 + 2a; 1 ^ + a; 2 2 )(6 2 + c 2 ^) , l ' 



c F( ) = ~2ft(7 2 + IpI 2 + Ui~ 2|p| ui + |p| ^ 2 - uiu 2 ) - C7a;| 

21 ''^ (7 2 + K-|p|) 2 )(7 2 + (^i-|p|) 2 + 2|p|^2-2a; 1 ^ + ^2 2 )(& 2 + c 2 ^)' l J 

At this point, we have all the elements required. There are two dominant contributions 
coming from the two products of internal propagators whose poles are nearly coincident, so 
in the following, we write twice the contribution corresponding to one of them (in particular, 
the one that has the singularity in p° = |p|, see Eqs. ([U]) and ([32]) below). Making use of 
the double spectral representation to do the Matsubara sum and the fact that the angular 
dependence of the integrand is trivial, one obtains 

tt«/ r.N 2 2Ar nir f ^P f°° duJ 2 1 f°° <M . . 

nliK, 0) = -g 2 N f 5^5 ab / — ^ / — — / —nfM 

6 J [Zn) J-oo Z7T lV q — UJ 2 J-oo LTX 

x \F x {uj x ,uj 2 ) + F 2 (u 1 ,u 2 ) + 8F 1 {ux,u 2 ) + SF 2 (u 1 ,u; 2 )]. (29) 



ImlT^V, 0) = -\g 2 N fq H^ ah J -^- jT ^-n' f (co) [F 2 (u, 0) + 5F 2 (u, 0)] , (30) 



Now, it easy to make the analytic continuation iv q — > q° + i0 + and to expand the imaginary 
part to lowest order in q°. The result is 

(27T) 

where we have used the property (5)Fi(ui,U2) + {S)F 2 (u>i + uj 2 ,uj 2 ) = 0. Thus, the explicit 
expressions of F 2 (u>, 0) and 5F 2 (u, 0) are required, 

F 2 {u, 0) = -p + {u, p) 2 = ^ -2 , (31) 



N c {(^ - \P\) 2 + Tf) 
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which, in the limit 7/ — > 0, can be replaced by F 2 (u), 0) = — 2tt5(uj — |p|)/7/ and 5F 2 (u!, 0) = 
2ti5(uj — |p|)/(-/V"cT/)- These substitutions give the final result to leading logarithmic order 

bnnSV.o) = -^g^ (1- ^) = -^4w^y (33) 



B. The gluonic contribution 

We only need to consider transverse gluons propagating in the side rails of the ladder, 
since the longitudinal ones do not propagate at hard momentum. We again modify the 
gluon propagator in order to include the effects of a thermal width 2^ g 

p(u, p) = i (- 7 ' 2 + -— f^— -) , (34) 

|P| \(w- |p|) 2 + T g (w+|p|) 2 + 7^y 

where we use the gluon damping rate at hard momentum 7 9 = a s N C T ln(l/g). The vertex 
equation which sums the ladders is shown diagramatically in Fig.|j. The normal vertex is 

rf abc (q,P, -P~q) = i/ o6c [2^V - S iJ P k - *V]. (35) 

We parametrize the effective vertex as 



rf c (q,P, -P~q)= if abc WVf k (p) + BS^p k + C8*p> + D^-). (36) 



We need the effective vertex only when the indices j, k are contracted with the transverse 
projectors 

A2w, -P~q)= Vj n m i r(q,P, ~P - q)KM = if*j?P? k {P)A- (37) 

After a bit of algebra the equation for the effective vertex reduces to 

p'A = 2p l + 2g 2 N c p 2 TY: I S^A(p + kyD\p + k + q)D\p + fc)-L (38) 

- J (27T) d \k.\ A 

where D t (p) is the propagator obtained from the spectral density of Eq. (|3~3|). Within the 
same previous approximations, one obtains 

11 



A2K + * q , ",; p) = Kf^Vl ( 1 + ^f^TT ^ 



,, !- ,.., ^ . , ( 39 ) 



With this result, following the same steps as earlier, it is easy to compute the imaginary 
part of the gluon polarization tensor. It reads 

n 2 T 2 N 

Imllf/V, 0) = - 9 -L^ q °^5 ab . (40) 

y 7s 

Finally, the sum of both contributions gives the correct value for the color conductivity 
a c = oJp/^fg, where u 2 = (g 2 T 2 /18) (2iV c + Nf) is the plasma frequency. 

III. SUMMARY 

In this work we have derived the color conductivity to the leading logarithmic order. 
The main point here has been the evaluation of the Kubo formula for the current-current 
correlator within the framework of thermal field theory in the imaginary time formalism. 
We have shown how to do the summation of the ladder vertex corrections, making use of 
some approximations based on the kinematical regime of the scattering, which corresponds 
to exchange of quasistatic transverse gluons of soft momentum, g 2 T < |k| < gT. We have 
also shown how to introduce a double spectral representation of three-point functions in 
order to do the Matsubara sums involved in this formalism. 

The possibility of computing other transport coefficients such as viscosity or electrical 
conductivity by a similar procedure deserves further consideration. This kind of trans- 
port properties follows from an integral over the differential cross section multiplied by the 
square of the momentum transfer k 2 . The corrections due to Debye screening and Landau 
damping arising at order gT are sufficient to render this integral infrared finite, scaling as 
g 4 log(k* /m c \) ~ g A log(l/g), where k* is a scale separating semihard and hard momentum 
transfers, restricted by gT <C k* <C T but otherwise arbitrary. So the log(l/g) comes from 
the sensivity to the momentum scale m e \ ~ gT, while the log(l/<?) in the color conductivity 
and the damping rate comes from a sensitivity to the magnetic scale g 2 T. This means that in 
order to attempt a similar computation of these transport coefficients, one must include into 

12 



the rungs of the ladder the longitudinal part of interaction as well as the Landau damping 
effects. The infrared sensitivity to the scale g 2 T should be compensated by a similar one 
coming from the hard damping rate. 
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+ 




FIG. 1. One-loop contribution to gluon polarization tensor. 





FIG. 2. Ladder diagrams contributing at leading order to gluon polarization tensor. The gluon 
line with the blob denotes the HTL resummed propagator. 





+ 




FIG. 3. Equation for the effective vertex which sums the quark ladder contributions. 






FIG. 4. Equation for the effective vertex which sums the gluon ladder contributions. 
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